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Abstract. We calculate the radiative properties for a linear dipole-coupled chain 
of qubits. Using the explicit energy eigenstates of the system, we find the radiation 
patterns for spontaneous transitions from the one-photon eigenstates to the ground 
state of the system. We show that depending on the excitation of a specific atom, the 
radiation tends to be focused either along or perpendicular to the chain. We conclude 
with a derivation of the total decay rate of the one-photon eigenstates, and find the 
interesting result that for systems where the photon wavenumber is not much larger 
than the interatomic spacing, up to 94% of the eigenstates are subradiant, that is, they 
decay significantly slower than a single atom in isolation. 
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1. Introduction 

There is considerable interest in the study of the entanglement properties and radiative 
dynamics of systems that are comprised of qubits, where the number of qubits is known 
exactly, with the aim of using such systems for quantum information processing. An 
understanding of the radiative properties of such a system is of great use since it 
provides information about decoherence and the radiation (transmission) channels in 
a multi-atom system, and how useful the spectrum of eigenstates are with respect to 
their persistence in time. Such information could be useful for the study of controlled, 
decoherence-free information transfer through the chain. 

The solution of a linear spin chain with a Heisenberg interaction was first given by 
Bethe [1, 2]. He analysed the properties of a linear spin chain in order to develop a theory 
of ferromagnetism. The system was assumed to have periodic boundary conditions, and 
he introduced a diagonalisation ansatz which allowed the eigenstates of the system to be 
obtained. The form of these eigenstates could be interpreted as "spin waves" travelling 
along the chain. Using this solution, he was able to obtain the partition function for the 
system. From the partition function, one can calculate many thermodynamic properties 
which describe the system macroscopically, such as the magnetisation [3]. Subsequently, 
a great deal of attention was paid by later authors to determining the energy eigenvalues, 
and in turn the partition function. The ground state and one-photon excitations received 
considerable attention, e.g. in [1, 2, 3, 4, 5]. 

Research followed in Bethe's direction, in considering various types of linear spin 
chains and investigating their macroscopic properties. The XY Heisenberg chain was 
solved exactly by Lieb et al. [4], for example. They also investigated correlations 
between spins, such as the long-range order. The anisotropic XYZ model was treated 
by [3, 5], who studied various macroscopic properties, including the free energy and 
susceptibility. The partition function and the energy eigenvalues of the XY chain were 
also presented and were in agreement with [4]. The ground state energy for the XYZ 
model was later obtained in the general case in [6]. 

A number of analytical studies relevant to quantum information processing have 
been performed on systems composed of A dipole-coupled atoms confined to fixed 
positions, with a particular interest in ring arrangements, i.e. a small number of atoms 
arranged in a ring at equidistant positions [7, 8, 9, 10]. Such forms require periodic 
boundary conditions and possess A-fold symmetry, thus making it easier to determine 
the energy eigenstates and energy levels. There is an introduced complication, however: 
in a line, one can consider only nearest-neighbour couplings; in a ring, additional 
couplings beyond the nearest neighbour must be considered. This clearly follows from 
the geometry of a ring as opposed to a line. Refs [7, 8, 9] analyse ring arrangements 
and develop a technique which can be applied to any A-atom ring system. It is 
also demonstrated that it is important to consider all interactions, not just nearest 
neighbours. As the solutions are not generalised for arbitrary A, a list of eigenstates 
are given for various values of A, but only up to the first three excitation levels. 
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The eigenstates corresponding to the elementary excitations of a ring and a ring 
with an atom at the centre are given in [10]. These states are used to study the associated 
decay rates. Energy eigenstates and eigenvalues have also been obtained for a linear spin 
chain of alternating atoms in [11]. Many papers have appeared on various aspects of 
entanglement in linear chains of qubits [12, 13, 14, 15, 16]. 

Recently, we have presented the analytic expressions for the energy eigenstates and 
energy eigenvalues for a linear chain of identical, two-level atoms (qubits) [17]. Such a 
system is mathematically equivalent to the Heisenberg XY model, which was analysed 
by Lieb et al. [4]. Although they diagonalized a more general Hamiltonian, we will 
demonstrate that our restriction to a specific interaction model, and not presenting 
solutions in terms of the Jordan- Wigner transform makes for a more perspicuous 
presentation in that the eigenstate basis is obtained using a new method, and the basis 
used is directly related to what can be measured, i.e. the excitation states of the atoms. 

A linear chain of atoms has been realized experimentally using ion trap techniques 
[18]. In these experiments the systems are controlled by individually interacting with the 
ions with focused laser light [19, 20]. However, the experimental form does not make 
use of the dipole-dipole interaction. Our approach, on the other hand, is to explore 
the dipole-dipole interaction in its potential to create multi-particle entanglement. 
Moreover, it allows one to selectively control atoms through the collective excitations of 
the system rather than individually addressing atoms with focused laser beams. To do 
this, however, we need to know how the system will interact with the electromagnetic 
field. 

In this paper, we use the state vectors given in our previous paper to explore 
the radiative properties of the system. We are specifically interested in the angular 
distribution of the radiation emitted by states that lie within the first excitation 
manifold. In addition, we will determine how many one-photon eigenstates are 
subradiant, i.e. exhibit reduced decay rates, and discuss its dependence on the number 
of atoms. In quantum information processing, it is ideal that we have a decoherence-free 
subspace. Although in general most of the eigenstates in the system presented here 
decay after a finite time, most of the one-photon states are subradiant, which is also 
beneficial. 

The paper is organized as follows. We first look at state properties of the chain in 
Sec. 2. In this section, we summarize the main result of the previous paper in which 
we derived the energy eigenstates and eigenvalues of the chain. However, we are now 
able to present the normalized states; the derivation of the normalisation constant is 
given in Sec. 2.1 on orthonormality relations. Next, we derive a parameter which gives 
a measure of the symmetry of a particular state, of which there are two possibilities: 
symmetric or antisymmetric with respect to reflection through the middle of the chain. 

In Sec. 3, we derive the transition probability per unit solid angle dT g /dVt in a 
transition from a M = 1 state to the ground state. The analysis concludes with a 
demonstration that the results for a two-atom system which are known in the literature 
can be recovered from our present work. We next explore in Sec. 4 the decay of the 
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Figure 1. As an example of our basis notation, the collective state shown in the figure 
for a 6 atom chain is represented by |3, 5). 

system when it is known which atom is excited, and find that the transition probability 
can be expressed as a simple superposition of eigenstate transition probabilities. In Sec. 
6, we look at the total decay rate of the one-photon eigenstates so that we can get an 
idea of the size of the subradiant subspace. 



2. General properties of the eigenstates 

Recently, we developed a new technique for finding the multi-quantum eigenstates 
in a linear chain of N identical, equally spaced and confined to fixed positions, two- 
level atoms each interacting only with its nearest neighbors through the dipole-dipole 
interaction [17]. The system we consider is governed by the Hamiltonian 

N N 

H = H + V = hu J2S? + ™ £ S+Sr, (1) 

i=l *>i=i 

where Hq is the interaction-free Hamiltonian and V is the dipole-dipole interaction 
between the atoms. 

In (1), ujq is the transition frequency of a two-level atom in isolation, Sf is the 
energy operator of the i-th atom, and Sf, are the raising and lowering operators for 
the i-th atom, respectively. The dipole-dipole interaction parameter fl = fi(ry) depends 
on the distance a between adjacent atoms in the linear chain and the polarization of 
the atomic dipole moments [x relatively to the interatomic axis. 

We can index the atoms in the chain by the numbers 1 to iV; collective excitations 
the system can then be represented using the ket \k\, . . . , kit), where the non-zero 
integers k\, . . . , ku denote the indices of the atoms which are in their excited state: 

\k l ,...,k M )=S+---S+JO), (2) 

where |0) denotes the ground state of the collective system, and is nondegenerate. A 
graphical example of this notation is given in figure 1. 

As we have shown in [17], the unnormalized eigenstates for the system are given by 

E C^\ kl ,...,k M ), (3) 

£:i<-"<£:m 
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with corresponding eigenvalues 

M 

4™ = M ^o + 2hQ £ cos (9i0, (4) 



=i 



where 



C 9v'-9m = e "-d M9ah0 ■ ■ ■ sin(g d k M 0, (5) 
and C, = tt/(N + 1). The integer 1 < M < N denotes the number of excited atoms in the 
system. The integers g±, . . . ,gM, where 1 < g\ < • • • < gu < N, label the eigenstates. 
In the present paper, we will work with the normalized form of the eigenstates: 

W-,J= £ c£::Jzl*i, ■■■,*">, (6) 

ki<---<k M 



where 



.\ M/2 



ki-k M — ( ^5.\ rik x -k M tn\ 
U 9i-9M I 7T / 91-9M ■ \' ' 

The derivation of the normalisation constant is given in Sec. 2.1 on orthonormality 
relations. Note the slight change in notation - unnormalized states, like those given 
in our previous paper, are specified with the M parameter enclosed by parentheses, 
whereas normalized states do not have the parentheses. 

2.1. Orthonormality relations 

It is clear that states with different excitation levels M are orthogonal. Here, we 
show that states with the same M but different gi are also orthogonal. Eigenvectors 
corresponding to different eigenvalues of a Hermitian matrix are automatically 
orthogonal. However, our some of our eigenstates are degenerate and hence it is useful 
to show that the eigenstates we use are orthogonal, no matter the multiplicity of their 
associated eigenvalues. As a byproduct of this proof, we also obtain the normalization 
factor for an arbitrary eigenstate as shown in (7). Consider the expression 

/ 7 // M ) L/,W \ \ " fik x -k Mr k x -k M (n\ 

ki<---<kM 

Although the C^ 1 .'.'.'^ coefficients are antisymmetric with respect to the k i: the 
product Cg/ L .'.'.'g' M Cgr--gM * s sy mm etric under such an interchange, for the two elements 
of the product reverse sign and hence the overall product does not change sign. This 
means that such a permutation will not affect the inner product of the two states; we 
introduce on the right-hand side of the above equation a sum over all permutations of 
the hi, which simply produces M\ copies of the inner product. To maintain equality, we 
divide this new sum by the factor Ml: 

cr(ki---kM) k\<--<kM 

The double sum in the above expression can be simplified. Permuting the hi merely 
implies that each ki can go beyond its range of k^i < ki < k i+1 . This means that we 
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can relax the restriction in the second sum that k\ < ■ ■ ■ < 1cm, and let each ki range 
from 1 to N. The only extra terms this new summation introduces are ones where two 
or more of the ki are equal. However, the Cgl'. '. '.g^ vanish in such cases, so equality is 
maintained in transforming the double summation into one where all the kj range from 
1 to N: 

fci,...,fc M =i 

Let us now consider the form of each term in this summation: 

C X : Z C ^ :: 9m = '"••••/""-/ sin 9aM sin g a kii • • • sin g),k M i sin g f k M ^ (11) 
In order the right-hand side of (10) to not vanish when summing over the ki, we 
must have g' a , = g a , . . . , g'f, = gj. This means that if the gi for one state are different 
when compared to another state, the inner product of those two states vanish. Only one 
individual g^ need be different in order to make the entire expression vanish. We now 
consider the case where the gi for the two states are the same, so that we may extract 
a normalisation constant. In the above expression, we see that each ki occurs twice in 
the form of two sine factors. It is clear that in the case where the & are identical for 
the two states considered, the only terms that contribute to a nonvanishing summation 
are those for which a particular ki is multiplied by two identical & in the two respective 
sine factors. A summation over each ki then produces a factor of (N + l)/2. Since there 
are M sine pairs, each nonvanishing summation evaluates to (N + 1) M /2 M . However, 
given one such nonvanishing sum, there will exist other nonvanishing sums of the same 
form but with the gi in a different order. This is due to the presence of the two Levi- 
Civita symbols. There are Ml permutations of the g^, so the summation will evaluate 
to M\(N + 1) M /2 M . Therefore 

<CJ«U = (^)\» - ^ (12) 

where it is assumed, as always, that g\ < • • • < gu- We therefore deduce that the 
normalized eigenfunctions are 




2.2. Symmetry properties 

One would expect that the energy eigenstates of a linear chain would be symmetric or 
antisymmetric when the line is reflected through the centre of the line. We now show 
that this is the case. To reflect a state through the midpoint of the line, we replace all 
the ki in the coefficients C^ 1 ".'*^ by k[ — N + 1 - fcj. We thus have 

k' ■■■k' i i 

Cgl- 9 M = tab -f sin g x k a i ■ ■ ■ sin g M k f £ 

= e ab ... f sm 9l (N + 1 - k a )i ■ ■ ■ sing M (N + 1 - k f % (14) 
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Next, we expand the sine terms as 

sin&(7V + 1 - kj)£ = sin g t (N + l)£cos&%£ - cos g^N + l)f sin&fcjf 
= sin g^n cos gikj^ — cos g^it sin gikj£ 
= (-l) 1+ ^smg i k j C. (15) 

Hence we have 

= e a6 ... / [(-l) 1 ^ 1 sin^U] ' ' ' l(-l) 1+9M sing M k f $ 
= (-l)M+9i+-+9« €ab ... f sin 9l k a ^--- sin g M k f ^ 

— (^ — \^ M +9i^ Y9M (jk\—k M (16) 

k' --k' 

We must remember that Cg\- g ^ is not in a canonical format, since the k\ are in the 
wrong order. To bring them in the correct order (i.e. complete reversal) requires 

(17) 

1=1 Z 

transpositions, hence there will be an extra factor of (— 1) M ( M_1 )/ 2 . 

If we think of an operator R which reflects a given eigenstate through the midpoint 
of the atom chain, we can say that it operates in the following way: 

*ICL> = (-i) m(m+1)/2+3i+ - +9m i« ) sm )' ( 18 ) 

where we have collected all the M terms together. Thus if (— i) m ( m + 1 )/ 2 +sh \-qm e q Ua j s 
+1, we say that the state is symmetric; if it equals —1, it is antisymmetric. This 
agrees with the result in the first paper, which says that for M — 1, gi-even states are 
antisymmetric and hence do not decay in the small sample model, since 1 + g\ is an odd 
number and hence (— 1) 1+91 = — 1. 

We would like to stress that the relationship between g and symmetry can be 
misleading for the case M — 1 due to an overloading of nomenclature. The words 
"symmetric" and "antisymmetric" are often associated with the words "even" and 
"odd," respectively, as is the case when discussing the parity of a function. One may be 
mislead into thinking that the symmetry or antisymmetry of an eigenstate is related to 
whether its g-value is even or odd. We must remember, however, that it is the even-g 
states that are antisymmetric, i.e. "odd states," and vice versa for odd-g states. 



3. Spontaneous radiation patterns for M = 1 to M = transitions 

In this section, we study the directional properties of the radiation emitted by the 
chain in the far-field zone (in the far-field zone, photons emitted from either end of 
the chain will differ in phase, but not direction). Directionality is an important factor 
as it determines the self-controlled directional properties of radiation that can be used 
to study the directional propagation of an excitation along the chain. When an atom 
is excited by an external field, it may spontaneously decay to its ground state. In 
the spontaneous decay of an isolated atom, the photon is usually lost and the atom is 
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ultimately left in its ground state. The radiation pattern produced by such an atom in 
free space is proportional to |^| 2 — |k-^| 2 . We will see, however, that the chain introduces 
an important modification to this pattern, which is due to the atom interacting with its 
nearest neighbours. 

Suppose that the system is initially prepared at t = in a single photon state 
(M = 1). We wish to calculate the photon emission probability per unit solid angle, 
dT/dfl, which gives the transition rate for a photon to be emitted in the direction (9, (ft) 
per unit solid angle and for the atom chain to go into the ground state (M = 0). The 
matrix element of V for the system in the M = 1 state to emit a photon of polarization 
e ks is given by 

1 /2 

where 

N 

fi = e ~ ik ' rj 57 + H.c. (20) 
i=i 

is the transition dipole moment operator and e ks is the polarisation vector of the emitted 
photon of wave vector k and polarisation s. The vector is the position of the i-th 
atom. We first calculate the transition dipole moment between the single-photon energy 
eigenstate and the ground state. Using (6), we find 

i=i 

N I 

., v , (kr 'sin.^. (21) 
3=1 ' 

By expressing the sine function in terms of complex exponentials, we find that there 
are two geometric summations generated by the sum over j. Using the formula for the 
summation of a geometric series, this results in the equation 

-ik-ri+ik-a • t 

where a = rj + i — r« is the displacement between two adjacent atoms. The transition 
probability per unit solid angle for a M = 1 state to emit a photon of polarization e ks 
is then given by 

dr *' s - ^M\v\o)\ 2 Vlj2 




f u 3 k \e ks ■ n\ 2 \ /sin 2 #A 1 - (-!)» cos[(7V + l)k • a] 



\ 87T 2 c 3 e h J\N+lJ (cosk-a-cos^) 2 
We now sum over the two polarizations s of the radiation field. By completeness, 
we have 

|e k i • n\ 2 + |e k2 • n\ 2 = \n\ 2 - \n ■ k| 2 , (24) 
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where k is the unit vector in the direction of the emitted radiation wave vector k. The 
total transition rate per unit solid angle in the far-field zone is then given by 

dTg ( uj 3 k {\n\ 2 - \k ■ n\ 2 )\ /sin 2 #A l-(-l)flcos[(iV + l)k-a] 




dfl \ 87i 2 c 3 e h J\N + 1J (cosk • a- cos#£) J 



(25) 



(26) 



\ 8n 2 c 3 eoh 

The factor in parentheses is the familiar spontaneous emission rate of one atom in 
isolation; the geometry of the chain and its multi-atom effects have been absorbed into 
the structure factor /(k, g, N): 

= l-(-l);co B [(iV + l)..al 

(cosk • a — cosgQ 2 \N + 1J 

The structure factor f(k,g,N) provides a measure of the collective behaviour of the 
chain, due to the specific arrangement of the atoms. We will focus on this object for the 
rest of the paper, to see what radiative effects are unique to a linear chain, and think 
of the structure factor as a "radiation pattern." It is important to remember, however, 
that the actual radiation pattern will also be modulated by the sine of the angle between 
the emitted photon's wave vector, and the atomic dipole moment. 

It can clearly be seen that the structure factor is symmetric about the centre of 
the chain, as evidenced by the fact that all the k • a terms are operated on by cosine 
functions. There is also the interesting factor 1 — (— l) 9 cos[(iV + l)k • a] present in 
the equation. In the small sample limit, this would become 1 — (— l) g , which would 
suggest that in such a limit all the g-even states would not decay, in agreement with 
our previous paper. Thus g being even is a sufficient condition for an eigenstate to be 
subradiant in the small sample limit. When g is odd, the state may be subradiant or 
superradiant depending on the particular case. We must keep in mind however that 
these classifications becomes more tenuous when phase differences between adjacent 
atoms become significant. 

One may notice that the structure factor as given above becomes undefined when 
cosk • a = cosg£. However, this ambiguity is only an apparent one. If we repeat the 
derivation of the detection rate formula, evaluating it before summing the geometric 
series, we obtain the following: 

N + 1 

f(k,9,N)\ cosk . a=cosg z = — g— , (28) 

a result which is consistent with taking the appropriate limit of (27) at the undefined 
point. A graph of the structure factor, (27), is plotted in figure 2 for a five-atom chain 
for a few g values. One can clearly see an overall decrease in the structure factor for 
higher values of g up to g = 3. Also, we note that the radiation is not isotropic, being 
either focused along the length of the chain, or perpendicular to it according to the 
value of g. 
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Figure 2. Plot of the structure factor f(k,g,N) for a five-atom chain against 
observation direction k in radians for different values of g. The cases plotted are 
for g = 1 (solid line), g — 2 (dashed line), and g = 3 (dotted line). 



The small sample behaviour of the structure factor is also interesting. Taking the 
limit k • a — > 0, we obtain the expression 

l-(-l) 9 /sin 2 c/A 



(l-cos^) 2 

l-(-l) 3 

7V + 1 



cot' \g^. 



(29) 



The form of (29) displays the radiation rate of closely spaced atoms. It is zero for even 
values of g, independent of the number of atoms. For odd values of g, the cot 2 term 
decreases as g increases from 1 to \N/2], and then increases from \N/2] to N. For a 
large number of atoms, the maximum value the structure factor can attain is 



7Y Z 



(30) 



while the lowest value for the structure factor is zero. 
3.1. Radiation from two and three atoms 

As an illustration of our general formula for the radiation pattern, consider a chain 
composed of iV = 2 atoms. In this case, we can recover the known form of the structure 
factor for two atoms as given in the literature [21, 22]. For N = 2, the structure factor 
(27) simplifies to: 



1 - (-l)flcos(3k-a) 
(cosk-a + U-l) 9 ) 2 



(31) 



Using the Chebyshev polynomial T 3 (x) = 4x 3 — 3x, we observe that cos 3a = 
4 cos 3 a — 3 cos a. This allows us to simplify the above expression (by dividing the 
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polynomial in cos k • a in the denominator into the polynomial in the numerator) to 

/(k,<7,2) = l-(-l)»cosk. a , (32) 

which reveals the radiation pattern of the symmetric (g = 1) and antisymmetric 
(g = 2) states of the two dipole-coupled atoms [21, 22]. Taking the small sample 
limit (k ■ a — > 0), we see that symmetric states decay twice as fast as an isolated atom, 
whilst antisymmetric states do not decay. 

For comparison, the structure factor for N = 3 is given by 

/(k, g, 3) = 1 + 2c\c\\\ - cosk • a - 2(cJ) 2 (-l)* cos 2k • a, (33) 

where c k g = \J2/(N + 1) sin gk£. In the small sample limit of ka 0, the structure 
factor is independent of the wave vector k and thus becomes a direct measure of the 
total collective decay rate. In particular, an eigenstate is superradiant or subradiant if 
/ > 1 or / < 1, respectively. We find that for the case of N — 3, in the small sample limit 
of (33), there are two subradiant states since f(g = 2) = and f(g = 3) = 3/2 — y/2, 
and one superradiant state since f(g — 1) = 3/2 + \[2. 

4. Radiation from one excited atom 

In the previous section we demonstrated the general radiative properties of a chain of 
collectively radiating qubits. We examine here the radiation pattern that arises from 
the transition \j) — > |0). The state \j ) corresponds to one atom in the chain being in 
the excited state whilst all the other atoms are in the ground state. In considering such 
a transition, we can see how the dipole-dipole coupling, or more pictorially, the presence 
of the other atoms in the chain, affect the radiation pattern of a single atom through 
the structure factor. The lowest order process that can occur is an initial state \j ) at 
t\ propagating to a time £3, emitting a photon, and then continuing to propagate to a 
time t 2 . For clarity, we set h — 1 in this section. 

The propagator for the nonradiating atom chain is given by 

K(2,l) = Y,e- lE ^- h) \^ g )(^\. (34) 
9 

The amplitude for the initial state to radiate a photon is given by 

(0\K V (2, > = i T dt 3 (0\K(2, 3)VK(3, ) 

= i^cj,(0\V\fl>l)2K5(u>-E g ), (35) 
9 

where we have set t 2 = —t\ = T/2 and approximated the integral over dts with a delta 
function. We now switch to a discrete electromagnetic field spectrum to convert the 
Dirac delta function to a Kronecker delta. This implies the replacement 

6(u-E g )^£-6 UtEa . (36) 
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To remain consistent, any integral over uj must be replaced by a discrete sum of the 
form 

^-^E- (37) 



I 



The expression on the right-hand side of (35) becomes 

<E4<°I^W> T ^.- ( 38 ) 

9 

We now square this expression to obtain a transition probability, and sum over all 
possible emitted photon frequencies u, weighted by the density of final states, to obtain 
the total transition probability per unit solid angle: 

-Jn T = T EE4^(^l^|0>(0|^|^)^ s ^pH. (39) 

w 9,9' 

Clearly, the sum over oo forces E g = E g > through the product of the two Kronecker 
deltas. Summing over g', we finally obtain 



= 27r5:(^) 2 |(^|^|0)| 2 p(^) 



an 

y 



9 



The radiation pattern from several different initial states for a seven-atom chain is shown 
in figure 3. Naively, we may expect two different radiation patterns in general: (a) the 
situation where the excited atom is at either end of the chain, and hence has only 
one interacting neighbor, and (b) the excited atom at any other position in the chain, 
surrounded by two interacting neighbors. However, as can be seen from the graph, this is 
not the case. This is probably due to the fact that the eigenstate excitation amplitudes 
(i.e. the c? g ) vary along the chain. Given the reflection symmetry of the chain, one 
instead classifies the patterns into \N/2] groups, since the j-th and (N + 1 — j)-ih 
atoms should radiate the same way by symmetry. Thus, the radiation pattern from the 
system depends in an essential manner on the initial excitation of the system. 

It is also interesting to see how the radiation pattern varies as a function of the 
number of atoms. This is shown in figure 4, where the initial state is |1) (one atom 
at the end of the chain excited, all others in the ground state), and the distance 
between adjacent atoms is held constant. We can see an interesting modification in the 
directionality of the radiation pattern with the number of atoms. Whilst the radiation 
pattern is spherically symmetric for two atoms, it becomes more concentrated along the 
axis of the chain. 



5. Diffraction properties of the structure factor 

In the preceding sections, we looked at how the structure factor shapes the radiation 
field for two different states of the chain. We now focus on general properties of the 
structure factor itself and how, by altering the number of atoms in the chain and the 
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Figure 3. Structure factor for a seven atom chain (ka = 1), for the initial states |1) 
(solid line), |2) (dashed line), |3) (dotted line), and |4) (dot-dash line). 
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Figure 4. Structure factor profile of the state |1) as a function of the number of atoms 
and angle of emission relative to the interatomic axis with ka = 1. 




distances between them, one may alter the directional properties of the radiation field 
emitted by spontaneous emission. 

From (27), we may expect to see effects in the radiation distribution similar to 
TV-slit diffraction. Using trigonometric identities, the structure factor can be expressed 
in the following form: 



f(Kg,N) 



sin 2 |(iV + l)(k.a-^) / sin 2 ^ 
2sin 2 |(k-a + ^)sin 2 i(k-a-^) 1^ + 1, 



(41) 
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Figure 5. Plot of the structure factor /(k, g = l,N) for a five-atom chain against 
observation direction k in radians for different values of ka. The cases plotted are 
ka — 0.1 (solid line), ka — 1 (dashed line), and ka — 8 (dotted line). 

We see the the structure factor contains the same modulating factor type as exists 
in the diffraction of iV + 1 slits [23] , namely 

sin 2 [|(iV + l)(k • a - gO]/ sin 2 [|(k ■ a - gft], (42) 

where the factor (k • a — g£) is analogous to the phase difference between two adjacent 
slits. 

We therefore expect that any diffraction maxima get narrower and more higher- 
order peaks are introduced as more atoms are placed in the chain. This can be seen if 
we plot the structure factor for a few different values of ka, as shown in figure 5. 

We can also analyse diffraction-like effects from single-atom radiation, shown in 
Figures 6, 7, and 8. In these graphs, there are several peaks which do not change as 
the number of atoms in the chain is increased. These are analogous to the the principle 
maxima for the diffraction of (N + 1) slits. We see that as more atoms are added to 
the chain, the positions of these maxima do not change, and that an increasing number 
of smaller maxima are added in between, consistent with our analogy of multiple-slit 
diffraction. All peaks do get narrower as the number of atoms is increased, which is 
consistent with our diffraction interpretation. 

Clearly, the structure factor profile is punctuated by peaks which become 
increasingly sharp as more atoms are added to the chain. In addition, reducing 
ka decreases the number of peaks. Therefore, to make the emitted radiation more 
directional along the interatomic axis, one simply reduces ka (or more physically, the 
distance between the atoms), and increases the number of atoms. When ka > 2n, 
one would expect repetitions in the angular profile of the structure factor. Thus, ka 
should be roughly 2tt. If it is significantly less than this, one starts to lose most of the 
interesting features of the structure factor, such as the sharp peaks shown in several 
previous figures. 




Figure 6. Structure factor of the state ||\ZV/2~|) for different values of N, ka = 0.1. 
The cases plotted are for N — 3 (solid line), N — 8 (dashed line), TV = 13 (dotted 
line), N = 18 (dot-dashed line). 




Figure 7. Structure factor of the state ||"7V/2]) for different values of N, ka = 2n. 
The cases plotted are for N = 3 (solid line) and N = 13 (dashed line). 



6. Total decay rate 

We can also find the total decay rate of a particular M — 1 eigenstate by integrating 
over all the possible directions of the wave vector k. This integration can be simplified 
by re-deriving the total differential decay rate (26) without step (22), the original 
purpose of which removes the summations in the expression, but makes integration 
harder. Following this modified procedure, we arrive at an alternate expression for the 
differential decay rate: 

dr, _ / ^(|^| 2 -|k-^| 2 ) \ » tk . a(p - q)p q , . 




Figure 8. Structure factor of the state ||~A72~|) for different values of N, ka — 6tt. 
The cases plotted are for N = 3 (solid line) and N — 13 (dashed line). 
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Figure 9. The decay rates of all eigenstates plotted against the number of atoms in 
the chain for the case /x • a = and k a = 1. Note that for this graph 7=1. 



Integrating over all directions gives the result 



N 



r s = 7 + 7 E F(k rij y 9 4, 



where 7 is the decay rate for a single atom in isolation and 



F(kr. 



u ) 




(A 



, On sin kv^n r 

J + [1 



kr 



3(A • iijf] 



cos kr; 



sin krij 



(kr i:j ) 2 (kr 



u > 



(44) 



(45) 



We have defined = r^/r^, = — rj, and = | 

We would like to get a sense of consequences of this equation. For example, some 
might expect that the collective interaction of the atoms through the vacuum field 
hastens decoherence of the system. Let us explore such a suggestion here. 



Radiative properties of a linear chain of coupled qubits 



17 



Figure 9 shows a plot of the decay rates for all M = 1 eigenstates as a function of 
the number of atoms. This plot is to be understood as follows: there are N eigenstates 
in the M = 1 subspace, where N is the number of atoms in the chain. These states are 
represented by where 1 < g < N; the decay rate of a particular eigenstate \ifj g ) is 
given by T g . Thus, for a chain with iV atoms, we will have N decay rates T g for the 
M — 1 subspace. Thus in the plot shown, there are N points associated with each value 
of N. 

We see that there is an interesting "banding" feature in the decay rates; the main 
features are an upper cluster of states or "band," a lower cluster of states, and a small 
number of states interspersed between the two clusters. The thickness of the "bands" 
varies with the orientation of the atom dipoles with respect to the chain axis. We 
can distinguish three groups according to the rate at which they decay: superradiant 
(r s > 7), natural (T g = 7) and subradiant (T g < 7). From these definitions, we 
may classify the upper and lower bands as superradiant and subradiant bands. The 
subradiant band is the most interesting feature in Figure 9. Since the eigenstates in this 
region decay slowly, such states have potential applications in regards to the storage of 
quantum information and photon trapping. 

To explore the size of the subradiant subspace, we consider the number of subradiant 
states against N, using the previously mentioned definition T g < 7 for such states. 
Figure 10 reveals that the number of subradiant states varies linearly with the number 
of atoms in an approximate fashion. For the few values of ka considered in the graph, 
the linearity of the relationship is maintained. Indeed, ka appears to have the general 
feature of only altering the gradient of the line of best fit, as can be seen in the graph. It 
is important to remember that we are dealing with a subspace of eigenstates, so Dicke's 
result [24] that half of the eigenstates are subradiant does not apply since that result is 
meant for the complete space of eigenstates. 

Now imagine lines of best fit which run through the roughly linear plots in Figure 
10. The gradient of such a line can be interpreted as a rough estimate of the fraction 
of the number of eigenstates which are subradiant; let's say that the gradient of the 
best fit line was 0.42 - then roughly 42% of the eigenstates are subradiant. This simple 
relationship is due to the fact that the line passes close to the origin, so we don't have 
to worry about a "non-zero y- intercept." Such a quantity is useful since it allows us to 
make general statements about eigenstate lifetime independent of the number of atoms 
in the chain, a consequence of the linear relationship. Obviously, such an approximation 
only works only when there are more than a few atoms in the chain. 

Let's now take a closer look to see how one can define a line of best fit. Plots 
like the ones shown in Figure 10 have the general form of straight lines with kinks in 
them to give a stepped appearance. We can compute the gradient by taking the average 
of the slope of the lines connecting adjacent points. This averaging procedure gives 
more weight to the straight lines that appear in the plots rather than the disjointed 
"steps" that connect them. Such an averaging procedure gives a good line of best fit for 
the ranges of ka and N considered here. The uncertainty in predicting the number of 
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Figure 10. Subradiant state count as a function of the number of atoms for the 
following cases: ka = 4.4 (marked by x's), ka = 1 (marked by -'s) and the small 
sample limit ka = (marked by +'s). We have set fi ■ a = 1. 

subradiant eigenstates from the gradient is of the order of a few states, but adequately 
serves the purpose of being a convenient semi-quantitative estimator for the size of the 
subradiant subspace for a linear chain. 

We now have a quantity which is independent of the number of atoms in the chain: 
the gradient of the subradiant count, or "subradiant fraction" (i.e. "the fraction of the 
total number of eigenstates which are subradiant") to be more concise. Let us have 
a look at how altering the only other major property of the system, ka, affects this 
quantity. 

A plot of the subradiant fraction (i.e. the "gradient of the best-fit line to the 
subradiant count," the "fraction of the total number of states which are subradiant") 
against ka is shown in figure 11. One can see that for small interatomic separations such 
that ka < 2, more than half of the eigenstates are subradiant, approaching 94% in the 
small sample limit. As ka increases beyond 2, we see a periodic pattern with interesting 
valley features where the number of subradiant states drops below 50 percent. 

The fact that most of the eigenstates are subradiant for ka < 2 is a beneficial 
property of this system. For example, to store quantum information, we need states that 
do not change in time for as long as possible. The ideal state would be one that does 
not decay. A collection of such states in which one can process quantum information 
is known as a decoherence-free subspace. Although this system does not have such a 
decoherence-free subspace, its subradiant subspace is very large, and such a collection 
of states is still useful for quantum information, given that the other approaches to 
quantum computing [25, 26, 27] do not rely on collective behaviour and hence are 
limited by the natural decay rate 7. 

Finally, we comment about the problem of individual addressing of the atoms in a 
possible realisation of single-qubit operations. In the case of direct individual excitation 
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Figure 11. Plot of averaged gradient or "subradiant fraction" against ka for /xa = 1. 

and addressing of the atoms, this may be realised in a microcavity situation, where atoms 
can be kept at large distances that are still comparable to the resonant wavelength [28] . 
Qubits in our system can be manipulated collectively, which clearly follows from the 
collective nature of the energy eigenstates. Regarding the problem of the encoding of 
quantum information given that targeting a specific qubit requires resolution of the order 
of a resonant wavelength, a potential fix for this problem is to introduce a Stark shift 
to increase the spacing between the two energy levels in each qubit such that diffraction 
no longer presents itself as a problem. This would also eliminate the dipole-dipole 
interaction, which is useful when direct single-qubit measurements are required. 

7. Conclusions 

We have studied the radiative properties of a chain of dipole-dipole interacting atoms. 
Assuming that only one photon is present in the system, we have calculated the 
component of the radiation emission pattern corresponding to collective interactions, 
and also the total decay rates to study the decoherence property of the chain. 

We were particularly interested in the directional distribution of the radiated 
field and the presence of the subradiant states that decay slower than a single atom 
in isolation. We have shown that the radiation tends to be focused either along or 
perpendicular to the chain, and that in some ways the emission pattern is analogous 
to multiple-slit diffraction. The total decay rates were also explored to see how many 
states are subradiant. We found that the fraction of eigenstates that are subradiant 
depends on the distance between the atoms. For small interatomic distances, we found 
that the majority of eigenstates are subradiant. 
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